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i—* ' Abstract 
5h ! 

Q_|. A new approach is introduced to study QCD amplitudes at high energy and com- 

\ paratively small momentum transfer. Novel cut diagrams, representing resummation of 

' Feynman diagrams, are used to simplify calculation and to avoid delicate cancellations 

CN ■ encountered in the usual approach. Explicit calculation to the 6th order is carried out 

to demonstrate the advantage of cut diagrams over Feynman diagrams. 

> : 
o\ . 

£3 : 1 Introduction 

^t- : 

The rapidity-gap events recently observed at HERA |I[ revived the community's interest 
in Regge-pole description of high energy scattering. The word 'Pomeron' seldom heard in 
recent years is once again found in the lexicon of experimentalists. It is perhaps then a 
good time to have a new look at the connection between perturbative QCD and Regge pole. 
Non-perturbative effect may be important for momentum transfer of the order of Aqcd or 
smaller, but we shall avoid it by going to a larger momentum transfer if necessary. 

It was proposed by Low and Nussinov 0] some years ago that the Pomeron may simply 
reflect a two-gluon exchange in QCD. Diagrammatic calculations have been carried out 
fH, f|, H, H, |7|, H in the leading log approximation to substantiate this proposal and to study 
other aspects of the high energy near forward scattering amplitudes. We shall follow the 
classical approach of implementing physical (s-channel) unitarity by summing up Feynman 
diagrams in the Feynman gauge [Q, [|]. However, other approaches are available. The use of 
physical gauge and dispersion relation || can sometimes make things more transparent. One 
can also emphasize on the complex angular momentum aspect by concentrating on t-channel 
towers and their interactions as guided by i-channel unitarity |J . We shall speak no further 
of these alternate approaches except to remark that our cut diagram method, to be discussed 
below, can be thought of an intermediate link between the classical Feynman diagram sum 
and the dispersion relation approaches. 

Due to the gauge and non-abelian nature of the theory, perturbative calculations are 
lengthy and complicated. Even by ignoring self energy and vertex corrections, as well as 
renormalization effects on the grounds that they will not alter the qualitative nature of high 
energy scattering we are trying to learn, even by ignoring quark pair productions as a first 
step, a complete QCD calculation can be carried out only up to the 6th order. In the much 
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simpler case of QED, an 8th order calculation || was reported to have taken sixteen months 
and two thousand pages to complete. In the case of QCD, it is much worse, a complete 
calculation of the 8th order is not available, though partial calculations have been carried 
out 0. 

What makes the calculation so complicated is the large number of diagrams that has to 
be tackled, and the inevitable cancellations between them. The computation of each diagram 
is fairly straight forward, though somewhat lengthy at higher orders. After the individual 
diagrams are calculated, care must also be exercised to add them up because of the presence 
of many delicate cancellations. Leading-log dependences on energy get subtracted away, 
complicated functions of momentum transfer also disappear. What emerges at the end is 
a product that is surprisingly simple. To the extent that it has been verified, high energy 
near forward scattering is described by multiple reggeized gluon exchange, supplemented by 
elementary gluon production off the reggeons and by s-channel unitarity [/J. 

This outcome has been verified up to the 6th order 0, [|, [5]. As mentioned above, many 
delicate cancellations take place in the sum to enable this simple picture to emerge. A 
complete calculation in the 8th order is not available, but if one assumes these cancellations 
which take place up to the 6th order also occur in the 8th and higher orders, then this 
reggeized picture has been verified up to the 10th order 0. This gives a strong support to 
the conjecture that it is true to all perturbative orders. 

While the conjecture is attractive, it is impossible to verify or refute without a new 
method to simplify the complicated calculations. We report in this article such a new 
method, in which high energy scatterings are computed via cut diagrams rather than the 
normal Feynman diagrams. These cut diagrams are not the Cutkosky cut diagrams. The 
sum of all cut diagrams here is equal to the sum of all Feynman diagrams, and not just 
their discontinuities. These cut diagrams can be regarded as a resummation of the Feynman 
diagrams in which many of the delicate cancellations encountered in the latter have been built 
in and explicitly avoided. Consequently, the simplicity of the final sum is revealed already 
in individual cut diagrams, and not masked by terms to be cancelled as is the case with 
Feynman diagrams. Besides, individual cut diagrams are easier to calculate than individual 
Feynman diagrams. 

In this paper we shall introduce the formalism of cut diagrams, as well as an explict 
calculation to the 6th order to demonstrate its effectiveness. We shall discuss quark-quark 
scattering throughout but allow the quarks to carry any SU(N) color charge. The result 
should be equally valid for gluon-gluon and other scatterings because high energy processes 
are insensitive to the spins of the colliding particles. We will study pure SU(N) QCD and 
ignore the production of quark pairs. If one is bothered by considering quark-quark scattering 
while ignoring quark productions, one can consider instead gluon-gluon scattering, and as 
mentioned above, the result would be the same. 

In a subsequent paper 0] we shall use the cut diagrams to discuss the reggeization 
conjecture in higher orders. 

In Sec. 2 the result of Feynman diagram calculations to 0(g 6 ) is reviewed. How these 
diagrams combine to give a sum much simpler than the individuals will be discussed in some 
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Figure 1: Quark-quark scattering in QCD up to the 6th order. The thick lines at the top 
and bottom of each diagram are the fermion lines, and the thin lines are gluon lines. 

detail. In Sec. 3, we begin to examine more carefully these cancellations by using sum rules. 
In Sec, 4, the method of cut diagrams is explained, and in Sec. 5, this method is applied to 
4th and 6th order calculations to demonstrate the savings effected by this method. 



2 Elastic scattering up to 0(g t 

Leading-log calculations in QCD, at high cm. energy y/s and comparatively small momen- 
tum transfer = \f&, are discussed, among others flUH, in the book of Cheng and Wu 
||. Citations to the original literature can also be located there. We review in this section 
the calculation of quark-quark scattering to the 6th order done in this standard way. 

The relevant Feynman diagrams are given in Fig. 1. The 2nd order diagram is labelled 
A, the 4th order diagrams labelled Bl and B2, and the 6th order diagrams CI to C21. This 
last labelling is identical to the ones used in Fig. 12.7 of Ref. P||Tl|. 

Under the interchange of the Mandelstam variables s = (pi +P2) 2 and u = {pi—p'2) 2 , the 
spacetime parts of these diagrams are either self conjugate, or one goes into another. For 
example, A^A, B2<->B1, C2<->C1, C15^C16, C17^C19, and C18<->C20, under s^u. 

Normalizing the Dirac spinors to uu — 1, and designating the fermion mass and gauge 
coupling constant as m and g, with (3 = g 2 /2ir, the T-matrix element T = (g 2 s/2m 2 )A4 
receives the following contributions from the individual diagrams to M. [3,4] ||12||: 

A = -ix-Gi 
Bl = -f3\n(se~™)l 2 -G 2 
B2 = +/3(lns)/ 2 -(G 2 + cG 1 ) 
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(2.1) 



Part of C21 has been combined with CI to form CI, and the remaining part of C21 has 
been combined with C2 to form C2. Note that the result for different diagrams is consistent 
with their s<-m character, under whose exchange the T amplitude responds with the swap 
se~ m <r->s. See App. A and B for a brief discussion of how these amplitudes are computed. 

The two-dimensional vector A in the transverse direction is the momentum transfer. The 
functions J n (A) and I n (A) are defined as follows: 



The functions I 2 , 13, J2, J3 are denoted respectively by I,Ii,K,I 2 in Ref. ||. The infrared 
divergences of these integrals can be regulated by a mass, either put in by hand or via 
the Higgs mechanism. This regulation discussed in the literature [3-9] does not affect the 
following discussions so we shall ignore it. There is however also an ultraviolet divergence 
in the integral defining J 2 (A), but it turns out that this function disappears in the sum of 
the sixth order diagrams so it causes no trouble either. 

The factors Gj in (2.1) are the color factors, with c = N/2 for SU(N) colors. As defined 
in 0, Gi, G2, G3, G4 are respectively the color factors for the diagrams A, Bl, CI, and C15. 
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They can be computed by an elegant graphical method H from the SU(N) commutation 
relation and identities 
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Hence the combination of color factors Gj given in (2.1) for the various diagrams remain valid 
whatever the color of the quark is, although Gj themselves would be different for different 
color of the quark. See App. A for a brief discussion on the computation of the color factors. 
The sum of all the terms in (2.1), from A to C20, is 



It is important to note the various cancellations taking place to make the sum (2.4) vastly 
simpler than the individual terms appearing in (2.1). For example, 

1. In the fourth order, the leading term proportional to Ins is cancelled out between Bl 
and B2 in the color amplitude proportional to G2, though not in Gi. 

2. In the sixth order, the leading Ins contributions to G4 from C15 to C20 also add 
up to zero. The expressions given in (2.1) are not accurate enough to deal with the 
subleading terms. The term in (2.4) proportional to G4 is obtained separately from 
the eikonal formula. 

3. As a result of these cancellations, the energy dependence and the SU(N) (or c) de- 
pendence of the Gi amplitude is (g 2 c\ns) m , and those of G2, G3, G4 are respectively 
g 2 (g 2 c\ns) m , g 2 (g 2 lns) m , g 4 (g 2 clns) m . These dependences can be summarized all at 
once by introducing a different notation for the color factors. We shall use the no- 
tation Fj j to denote a color factor with i parallel vertical lines connecting the two 
fermions, and j parallel horizontal lines joining any two of the vertical gluon lines. 
We shall also write F i simply as Fj. The relations with the color factors Gj are 
Gi = F 1; G 2 = F 2 , G 3 = F 2| i, and G 4 = F 3 . The g, c and Ins dependences of Fy in 
(2.4) are then given by g 2 ( J_1 ) (g 2 c In s) m c~- ? for a diagram of order 2(m + i). We shall 
refer to such dependences as Regge-like, for the Reggeization of the scattering ampli- 
tude to be discussed later relies critically on this feature of the scattering amplitude. 
Note from (2.1) that contributions from individual diagrams are not Regge-like. Only 
the sum is. 
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a(A) = (3cA 2 I 2 (A) . 
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4. Simplification in transverse-momentum dependences also occur in the sum. The simple 
integrals I n survive, but the complicated integral J 3 and the divergent integral J 2 do 
not appear in the sum. This cancellation is highly nontrivial because both of them 
contribute different amounts to different color amplitudes. More specifically, 

5. The function </ 3 (A) appears in all the color amplitudes Gi, G 2 , G 3 and G 4 in diagrams 
C15 to C20. Those in G 2 , G 3 , G 4 actually get cancelled out in the sum, but its presence 
in the Gi amplitude survives. However, since this term is of order g 6 In s, it is negligible 
compared to terms of order g 6 \n s appearing in the Gi amplitudes of C2 and Cll to 
C14, it can be ignored in the leading- log result displayed in (2.4). 

6. J 2 (A) appears in the color amplitudes Gi,G 2 and G 3 in individual diagrams CI to 
C14 and all these appearances get cancelled out. 

As a result of these cancellations, M. acquires a very simple interpretaion in terms of reggeized 
gluon exchanges . These exchanges are constructed in such a way to ensure s-channel 
unitarity [3,7]. 

Let us denote the reggeon propagator by 

-Ri(A, s) = — exp(-a<(A)lns) . (2.6) 

This reduces to the (transverse part of the) ordinary propagator ii(A) = A -2 for small 
g 2 c\ias. Similarly, let us denote the reggeized version of I n (A) by 

R n (A,s)= J (fl|^^(^^)) ■(2n) 2 6 2 (±q i± -A > j , (2.7) 

indicating the exchange of n reggeons. Then to order g 6 in T, we can write 

a 2 

M = - J R 1 (A,s)-F 1 + i|-[ J R 2 (A, S )-F 2 + J R 2jl (A, S )-F 2)1 ] 

+ |^ 3 (A, S )-F 3 . (2.8) 

In other words, the Fi, F 2 , F 2j i, and F 3 components looked precisely like diagrams A, Bl, CI, 
and C15 respectively, but with the vertical gluons replaced by their reggeized version whose 
propagators are given in (2.6), and with all longitudinal-momentum integrations omitted. To 
interpret it this way for -R 2j iF 2j i we need to know the Lipatov-Dickinson vertex || describing 
how elementary gluons are produced and absorbed from the reggeized gluons. 

This remarkable simplicity and regularity led to the conjecture [7| that the reggeized 
formula (2.8), suitably generalized, is the correct high energy limit to all perturbative orders. 
This conjecture is very difficult to verify on account of the shear complexity in higher order 
calculations. For QCD in the 8th order it is simply not manageable without simplifying 
assumptions. If one assumes all cancellations occured up to 0(g 6 ) will also occur in higher 
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orders, the final result can be extracted from a relatively small set of diagrams, then it is 
reported that this reggeization conjecture is true to 8th and 10th orders H. Even so these 
heroic calculations are so lengthy and complicated that to our knowledge the full details 
have never been published. 

3 High energy kinematics 

We shall discuss in the next two sections a method of using cut diagrams to sum up the 
Feynman diagrams, a method in which most of the cancellations discussed in the last section 
are automatically built in. This simplification shortens the computations and makes it 
possible to study higher order diagrams. In this paper we shall discuss how it simplifies the 
calculations up to 0(g 6 ). In a subsequent paper ]!(]] we shall discuss how it helps to verify 
part of the reggeization conjecture to all orders. To prepare for the grounds for both, we 
discuss here the relevant kinematical features of high energy scattering which enables this 
new method to work. 

We will assume the colliding beams in their cm. system to be directed along the z 
direction. In lightcone coordinates, p± = p° ± p 3 , the components of a four- vector are 
labelled in the order p^ = {p+,P-,p±), with the two-dimensional vector p± lying in the 
transverse x-y plane. In this notation, the incoming fermion momenta are p\ = (a/s, 0, 0) 
and P2 = (0, v^i, 0), in which their mass m has been neglected. The outgoing fermion 
momenta are approximately given by p[ = (\/s, 0, A) and p' 2 = (0, y/s, —A). See Fig. 1. 

Suppose n gluons are hooked up to the upper fermion line as shown in Fig. 2. The initial 
and final fermions are on-shell but the gluons can be off-shell, though with an amount of en- 
ergy far less than ^/i. At high energy, the numerator of the propagator can be approximated 
by 

7P = 2m u x {pi)ux{pi) (3.1) 

A 

provided the Dirac spinors are normalized to u\(p)uy(p) = 5\y. With that, the dominant 
current ma(Pi)7°ma'(pi) at high energy is just its translational part 5x\'Pi/m. This shows 
that the spin content at high energy is unimportant. All that it does is to enforce helicity 
conservation of the fermion, and to produce a factor 2p± at each vertex together with an over- 
all normalization factor of l/2m. For most of the discussions below we shall ignore this QED 
factor and concentrate on the contribution from the denominators, i.e., the corresponding 
scalar theory. 

The denominator of the ith inverse propagator is 

i i 

(Pi + Qi) 2 -m 2 + it — s(^2 Xj + it) (3.2) 
i=i " i=i 
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Figure 2: A quark — > quark + n gluon tree diagram, 
where Xi = q^/ yfs. The scalar amplitude in Fig. 2 is given by s~ n times 

(n \ n— 1 i 
E g i n V t T . , v • (3-3) 
3=1 / i=i 2^j=i -fj -r it 

Note that a momentum conservation 5-function for the negative components (together with 
an explicit factor —2iri) has been incorporated. In eq. (3.3) the ordering of the gluon lines 
from left to right is [123 ■ ■ - n]. If they are ordered differently, say [v\V2 ■ ■ -v n ] = V, then the 
corresponding amplitude is 

(n \ n—1 j 

■ n ^ — ~~ • (3- 4 ) 
3=1 / i=l ^3=1 X f j + 26 



4 Sum Rules 

It is possible to compute sums of Feynman diagrams without much of the delicate cancella- 
tions discussed in Sec. 2, by using the cut diagrams which we shall describe in this section 
and the next. The derivation of cut diagrams relies on two exact combinatorial formulas for 
the quantity a[V] in (3.4): the factorization formula and the multiple commutator formula 
derived in Ref. [jl~4j| . We shall discuss the former in this section, and the latter in the next 
section. 

Consider an ordering of n« gluon lines: \vi\Va. • • -Vi ni } = V{. We shall use the notation 
{V\] V%\ • • • ; V m } to denote the set of a/Zorderings of the M = Y%Li n i gluon lines, providedthe 
relative orderings of lines within each Vi are maintained. The number of orderings in this set 
is given by the multinominal coefficient Ml/ HJLi n y- For example, if V\ = [135], V% = [24], 
then {Vi,V 2 } = {135; 24} consists of the 5!/3!2! = 10 orderings [13524], [13254], [13245], 
[12354], [12345], [12435], [21354], [21345], [21435], and [24135]. 

We shall use the notation 

a{V l -V 2 -----V m }= Y, a l V ] 

Ve{Vv,V 2 ;-;V m } 
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to denote the sum of all amplitudes for the gluon orderings in the set. The factorization 
formula MM then states that 



a{V 1 ;V 2 ;---,V m } = l[a[V i \ (4.2) 



i=i 



In particular, if each set Vi = [vj\ consists of only one gluon line labelled by Vi, then 
{Yi'i V2] " " " 5 Vm} is the set of all orderings of the m gluon lines. In that case the factor- 
ization formula reduces to the well-known eikonal formula |JJ, |T5 |. Other special cases of this 
formula have also been discovered before |7|, |T6" |. 



It is useful to adopt an alternative notation for the right hand side of (4.2) to denote 
a>[Vi\ simply as a[Vi| V2I • • • \ V n ]. This notation is suggestive because the vertical bar 
can be interpreted graphically as a cut in the fermion propagator between the last gluon 
line of Vi and the first gluon line of V i+ i. For a cut propagator, instead of the usual factor 
(J2)=i x vj +2e)~\ we have — 27ri5(£)*-_ 1 x v 3 )- This notation is also convenient because it makes 
(4.2) deceptively simple. It now reads a{Vi \ V 2 : • • ■ ; V m } — a[Vi[X^| • • • \V m ); we simply have 
to change the semicolons to vertical bars. 

Cut propagators are not limited to tree diagrams like Fig. 2. The off shell gluons can be 
connected to other diagrams to form a composite diagram that inherits the original cuts. 
The cut diagrams so formed are similar to but different from the Cutkosky cut diagrams. 
Similar because we have the same factors for the cut propagators. Different because the cuts 
here occur only on fermion lines whereas in a Cutkosky diagram they can occur on any line. 
Moreover, via (4.2), our cut diagram represents a sum of M\/ Yl™ =l rij\ (uncut) Feynman 
diagrams, with their real and imaginary parts fully included, which is unlike the Cutkosky 
diagrams in which only the imaginary part or the discontiuity is represented. 

It is clear from (4.2) that the factorization formula can be thought of as a sum rule, to 
represent sums of Feynman diagrams as cut diagrams. As will be discussed in Appendix 
B, a cut diagram is easier to compute than an uncut diagram. In this way not only it is 
unnecessary to compute the individual diagrams first, the cut diagram representing the sum 
is actually easier to compute than just one single Feynman diagram. 

We shall now apply the factorization formula to compute sums of amplitudes quoted in 
Sec. 2. We shall see that the Ins factor and the Jj functions that get cancelled out in the 
sum never once appear in the cut diagrams. 

In what follows we shall use the notation (Bl) to denote the spacetime part of Bl 
(without the color factor G2). Similar notation will be used for the spacetime part of other 
diagrams as well. 

In the 4th order, the G2 coefficient in M. is given by (Bl) + (B2). Using the eikonal 
formula a[12]+a[21] = a[l\2] on the upper fermion line (see Fig. 3), (Bl) + (B2) is reduced to 
a cut diagram that can be easily calculated (see Appendix B for all the calculations), yielding 
the correct result g 2 il2/2 given in (2.1) and (2.4). Similary, the sum rule a{l; 2; 3} = a[l|2|3] 
applied to the six horizontal ladder diagrams in the 6th order, yields a cut diagram (see 
Fig. 4) which gives the sum of (C15) to (C20) to be g A Is/6. This is the correct coefficient 
of G4 given in (2.4). Note that neither Ins nor J 3 ever appears. 
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12 12 12 

Figure 3: An illustration of summing Feynman diagrams to obtain cut diagrams. 
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Figure 4: Another sum rule. 



321 312 231 31,2 




123 123 123 123 



Figure 5: Yet another sum rule. 



3 2 1 2 1 3 1 3 2 2 1,3 2 3,1 1 3,2 




123 123 123 123 123 123 



Figure 6: An illustration of how sum of Feynman diagrams can be turned into sum of cut 
diagrams. 
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Let us now look at the cancellation of Ja(A) in the G 2 and G 3 color factors, as discussed 
in point (5) of Sec. 2. According to (2.1), the coefficient of — G3 from diagrams C15 to C20 
is given by (C16) + (C19) + (C20). Using a{31;2} = a[31|2] on the upper fermion line as 
in Fig. 5, this sum is given by a cut diagram which can be evaluated to be — g 4 i Ins 1^/271. 
Again J3 does not appear in the cut diagram. Similarly, the coefficient of CG2 from diagrams 
C15 to C20 is given by (2.1) to be 3(C16) + (C17) + (C18) + 2(C19) + 2(C20). This is 
equal to twice the previous sum, plus (C16) + (C17) + (C18). To use the sum rule, this 
requires an expression (see Fig. 6) for a[321] + a[213] + a[132]. It is easy to see that this 
is equal to a{21;3} — a{23; 1} + a{13;2}, and by the factorization formula, also equal to 
a[21|3] — a[23|l] + a[13|2]. There may seem to be very little gained by replacing the sum 
of three terms by the sum of another three terms, but this is not so. We are replacing 
the sum of three uncut diagrams by the sum of three cut diagrams. The computation of 
cut diagrams is much easier than the computation of uncut diagrams. For one thing the 
complicated function J 3 that appears on all the three uncut diagrams but disappears from 
their sum never appears in any of the three cut diagrams. Evaluating the cut diagrams, 
the coefficient of cG 2 from C15 to C20 is — g 4 ic\n sl 3 /2ii, agreeing with the answer given in 
(2.4). 



5 Cut Diagrams 

The method introduced in the last section has a serious shortcoming. It does not tell us 
which cut diagrams to compute without detailed considerations of the kind carried out 
there. In this section we discuss a remedy for this shortcoming with the help of the multiple 
commutator formula [12]. 

So far we have been treating Fig. 2 mostly as a QED or a scalar amplitude. For QCD 
the non-abelian color matrices t a have to be incorporated. Instead of (3.4) the amplitude is 
now 

A[v!V 2 ■■■v n }= a[viv 2 • --Vnjtlv^ ■■■v n \ = a[V]t[V] = A[V] , (5.1) 

where t[V] = t vl t V2 ■ ■ -t Vn . What we want is a formula for the sum of the n\ permuted gluon 
orderings, A = Y,ves n A\V\. For QED, where we can take all t a = 1, this is simply the 
eikonal formula, so what we need is the non-abelian generalization of it. This is the multiple 
commutator formula 

a= E *\y\t\y]= £ a[v c ]t[K]- (5.2) 

ves n v&Su 

It expresses the sum of a[V]t[V] in terms of sums over the corresponding cut amplitude 
afVj^Vj]. Compared to the eikonal formula this looks complicated; instead of a single term 
on the right hand side we have now a sum over n\ terms. The complication is inevitable 
because we are attempting to sum up amplitudes for every color. However, we shall see 
that many of these terms are actually zero, and moreover, the cut diagrams on the right 
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are considerably simpler to evaluate than the uncut diagrams on the left. Again delicate 
cancellations will largely be avoided as before. 

What remains to be described is what the cut diagram V c that corresponds to the 
Feynman diagram V is, as well as what the amplitudes a[V c ] and t[V^] are. Given a 
V = [v iv 2 • • • v n ] , start from the rightmost number v n and proceed leftward until one comes 
to the first number less than v n . Put a cut just to the right of this number. Then start 
from this number and proceed leftward again until one comes to the first number that is less 
than this number, and another cut is put just to the right of this new minimum number. 
Continue this way until the end and we have constructed the cut diagram V c . For example, 
for n = 2, the 2 cut diagrams are [12] c = [1|2] and [21] c = [21]. For n — 3, the six cut 
diagrams are [123] c = [1|2|3], [213] c = [21|3], [312] c = [31|2], [132] c = [1|32], [231] c = [231], 
and [321] c = [321]. 

To each cut diagram we associate a spacetime cut amplitude a[V c ] as described in the 
last section. Namely, it is given by (3.4) except the propagator at a cut is replaced by 
-2TriS(J2jX Vj ). 

The complementary diagram V' c of a cut diagram V c is obtained as follows. If a cut 
appears between two numbers in V c , then there will be no cut between the same two numbers 
in V^, and vice versa. For n = 2, the complementary cut diagrams are [1|2]' = [12] and 
[21]' = [2 1 1] . For n = 3, the complementary cut diagrams are [1|2|3]' = [123], [21|3]' = [2|13], 
[31|2]' = [3|12], [1|32]' = [13|2], [231]' = [2|3|1], and [321]' = [3|2|1]. 

When no cut appears in V' c the color factor t\V[\ is simply t[v\V2 • • • v n ] = t Vl t v , 2 • • • t Vn . If 
a cut appears between vi and Vi + i, then the product t Vi t Vi+1 is replaced by their commutators 
[tvi,t Vi+1 \- If two or more consecutive cuts appears, then the corresponding product of t's 
is replaced by multiple commutators. For example, t[2|13] = [t2,ti\ts, t[2|3|l] = [t 2 , [t3,£i]], 
andt[4|3|2|15] = [t 4 ,[* 3 ,[*2,*ip5- 

6 Cut amplitudes to 0(g 6 ) 

We shall compute the cut diagrams to 0(g 6 ) to demonstrate the simplifications obtained 
therefrom. 

Using (5.2) on the upper fermion, we obtain a set of cut diagrams as shown in Fig. 7. 
This set is not unique as others can be obtained from a different labelling of the gluon lines. 

In the horizontal ladder diagrams Bl, B2, C15-C20, we choose the planar diagrams Bl 
and C15 to be the ones whose upper lines are completely cut. This fixes the labelling [12] 
for Bl and [123] for C15 as shown, and it in turn determines how diagrams B2 and C16 to 
C20 are to be cut. Using computational methods discussed in App. A and B, these diagrams 
yield 

Bl c = +^g 2 il2-G 2 
B2 C = +/3(lns)7 2 -cGi 
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Figure 7: cut diagrams up to the 6th order. 



C15 c 

C16 c 
C17 c 
C19 c 
C20 r 



6 

-/5 2 (lns)2J 3 -c 2 G 1 

C18 c = 



-//3i(m S )/ 3 -(cG 2 -G 3 ) 



(6.1) 



Similarly, we choose to cut the line of the planar diagram CI to obtain 



Cl c =-g 2 (3i{\ns) 
C2 C = -/3 2 (ln 2 s) 



J2I2 



±A 2 /f - J 2 / 2 



•G 3 

C 2 Gi 



(6.2) 



The twelve diagrams C3 to C20 can be divided into four groups of three, each giving 
identical contributions, so it is necessary to consider only one of these groups. The group of 
C3, C7, C12 have been chosen for that purpose. There is a symmetry between gluon lines 
2 and 3 so we may double this group and consider it as a sum of six Feynman diagrams. 
By applying the multiple commutator formula, the six cut diagrams shown in Fig. 7 are 
obtained. Their values are 

C3 C = ^(C3 / c + C3 / c , ) = -^(lns)ij 2 / 2 -G 3 



C7 C = 1 -(C7' C + C7:) = 
C12 c = l -{Cl2> c + Cl^)=p\\n 2 s)- A J 2 I 2 .{-c>G 1 ) 



(6.3) 



The expressions in (6.1)-(6.3) should be compared with those of the uncut diagrams, eq. (2.1). 
It should also be compared with the sum found in eq. (2.4). Several points can be noted 
from these comparisons: 
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1. Ins factors that get cancelled in the sum of the Feynman amplitude (see points (1) to 
(3) in Sec. 2) never even appear in (6.1)-(6.3). Cancellation of this kind is automatically 
built into the cut diagram formalism. 

2. The tranverse function J 3 appears only in C16 c in eq. (6.1). This expression survives 
the sum but can be ignored compared to the contribution from C2 C . In other words, 
as opposed to the Feynman amplitude (2.1) where J3 appears in many places but most 
of them are cancelled out at the end (see points (4) and (5) in Sec. 2), in the cut 
amplitude J 3 does not appear except when it survives the sum. 

3. The cut amplitude is not as successful in cancelling the transverse function J 2 (point (6) 
of Sec. 2), although there is still an imporvement here over (2.1) in that J 2 appears in 
fewer places. In fact, it appears in C3 C to C20 c only when absolutely needed to cancel 
its previous appearance in Cl c and C2 C . In order for J 2 to disappear completely it 
is necessary to combine diagrams with triple and four gluon vertices together using 
the Lipatov-Dickinson vertex |J. The technique of cut diagrams by itself, which deals 
mainly with the fermion lines, is not sufficient for that purpose. 

4. Other than the J 2 complication mentioned above, the summands of the final answer 
(2.4) appear directly in the cut amplitudes. In that sense the cut amplitudes are as 
economical and as simple as they can ever be. In particular, the Regge-like feature 
mentioned in point (3) of Sec. 2 is present already in individual cut diagrams. 
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A Color factors 

The method to compute the color factors graphically || is briefly reviewed here. It is suitable 
for both cut and uncut diagrams. 

The basic tool is the graphical identities depicted in Figs. 8 and 9, which are nothing 
but the commutation relation and identities in (2.3). The commutation relation is valid for 
all representations of the color matrix t a , hence Fig. 8 remains true when the quark line is 
replaced by a gluon line. 

As an example, the computation of the color factor for the uncut diagram C19 in eq. (2.1) 
is carried out in Fig. 10, and the computation of the color factor for the cut diagram C19 c 
in eq. (6.1) is carried out in Fig. 11. For the latter, we have to remember to use the 
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Figure 8: Graphical represenation of the commutation relation of color matrices. Thick 
lines are fermions and thin lines are gluons. The color factor at a gluon-fermion vertex is t a 
of eq. (2.3), and the color factor for a triple gluon vertex whose color indices a,b,c in the 
diagram are in clockwise order is if a bc- A cut represents a commutator between two color 
matrices. 




Figure 9: Graphical represenation of the last two identities in equation (2.3). c = N/2 for 
SU(N) color. 




= G 4 + 2 c G 2 - G 3 

Figure 10: A sample calculation of the color factor of a Feynman diagram. 
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Figure 11: A sample calculation of the color factor of a cut diagram. 






(a) 



(b) 



(c) 



Figure 12: Flow diagrams. 



complementary cut diagram for the color factor, and not the cut diagram itself shown in 



The color factors in eqs. (2.1), (6.1)-(6.3) can all be obtained in this way. 

B Spacetime amplitudes 

High energy limit of a Feynman diagram can be computed either in the Feynman-parameter 
representation^], or directly in momentum space using lightcone coordinates^, |], [5]. We 
follow the latter approach and review briefly the main points. For a more detailed discussion 
see Ref. ||. It serves well to remember at this point that the bottom fermion carries mostly 
the ' — ' momentum and the top fermion line carries mostly the '+' momentum. 

The procedure to follow for the computation is simple: (i) carry out the '+' momentum 
integration using residue calculus; (ii) then carry out the ' — ' momentum integration to 
produce the Ins dependence. The transverse momentum integrations are never explicitly 
carried out, which is why the answers are all expressed in terms of functions like I n (A) and 
J n (A) given in eq. (2.2). 

Step (i) requires a knowledge of the location of poles for each propagator, which in turn 
depends on the choice of loop momenta k a . Suppose q = P + J2a=i c °^a is the momentum 
flowing through a propagator, with P being some combination of external momenta and c a 
some integer coefficients. The propagator D~ x = (q 2 — m 2 + ie)^ 1 = {q^q+ — q± — m 2 + ie)^ 1 
has a pole in the integration variable k a+ if c a ^ 0. This pole is located in the upper half 
plane if c a g_ is negative, and in the lower half plane if it is positive. A way to keep track 
of the sign is to draw a flow diagram for the '— ' momentum, with arrows indicating the 



Fig. 7. 
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direction of the ' — ' flow. See Fig. 12 for some illustrative examples. The arrows around a 
loop pointing one way (clockwise or counter clockwise) have their poles in one half plane, 
and those pointing the other way have poles in the opposite half plane. The arrows in the 
flow diagrams must obey momentum conservation, and one may also assume they do not 
go around in closed loops. Otherwise the poles will all be in the same half plane, and the 
integral is zero if the contour is closed in the opposite half plane. With these constraints we 
see in Fig. 12 that a one-loop diagram allows only one flow path, but for two and more loops 
there is bound to be more than one flow diagrams because these constraints simply cannot 
fix the direction of the ' — ' flow on a boundary line of two loops. 

Each flow diagram corresponds to a range of k a - variables. By definition, the ' — ' variables 
along the direction of the flow are always non-negative. 

We shall always close integration contours in the lower half planes, and indicate the poles 
so enclosed by a cross (x) in the flow diagram. For a scalar diagram, the T-matrix element 
is equal to the product of / propagators D^ 1 = (q 2 — m 2 + ie)" 1 , integrated over the £ 
loop momenta d 4 k a = ^dk + dk_d 2 k±, with an extra numerical factor — [i/(27r) 4 ] £ . Each '+' 
integration produces a factor —2ni, the T-matrix is equal to 

T=-Y [ Dk ± (f]—)^- , (B.l) 

where x a = k a _/^/s is the scaled '— ' momenta, 

M -Sl? (b - 2) 

is the measure for transverse momentum integration, and D{ is either the propagator eval- 
uated at the x poles or the residue of the x pole divided by y^s. Each summand in (Bl) 
corresponds to a flow diagram with a pole taken from the lower half plane (an x pole) of each 
loop. A flow diagram may have more than one set of x poles, in which case the sum is taken 
over all possible sets. For example, the one-loop diagram Bl in Fig. 1 has only one flow path 
(Fig. 12(a)) and one set of x poles . The two-loop diagram C17 of Fig. 1 has two flow paths, 
Figs. 12(b) and 12(c). The first flow diagram has two sets of x poles, and the latter flow 
diagram has one set of x poles. We shall see later that in the leading log approximation, 
we may discard the q + dependence on the upper fermion line and the accompanied poles, in 
which case 12(b) is also left with one set of x poles. 

Equation (Bl) is also valid for cut diagrams, provided the propagator factor D^ 1 = 
(q 2 — m 2 + ie)" 1 of a cut line is replaced by —27ri5(q 2 — m 2 ). 

Let = % — 1, 2, • • • , I. Every z,i is a linear combination of x a , and in the case of 

a propagator along the bottom fermion line, x = 1 is also involved. The sign of each z% is 
fixed by the direction of arrows in the flow diagram. If there are x cu t lines, then the last x 
x a will be chosen to be equal to \zi\ of these cut lines. The / indices i will now be divided into 
three sets, a, b for indices from 1 to £ labelling the internal lines with an x pole, u for indices 
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labelling propagators on the top fermion line, and m for the rest. Then q a+ = cfci_/y/sz a , 
and q m +, q u+ can be expressed as linear combinations of q a+ , e.g., 



Qm+ — C ma q a+ — Z m c mbQb±/ Z b (B-3) 
a b 

Within the leading log approximation, one has 

D a = z a 

D u ~ sz u {uncut) , D u = —27ii5(sz u ) (cut) 

D m = C mb — ql ± - q m ± + ie ■ ( BA ) 

b Zb 

The ' — ' momentum flows mainly along the bottom fermion line, with very little seeping 
out to avoid a substantial mixing with the '+' momentum coming from the top fermion 
line, for a finite mixture of these two at a propagator would make it proportional to s and 
therefore negligible. This means the dominant contribution to T comes from regions where 
all x a are small. Since x a = for a = £ — x+ 1, • • • , £, on account of the (^-functions of the cut 
lines, the region of '— ' integration can be roughly divided into regions where the remaining 
— £ ~ X ^a's are strongly ordered, and regions where two or more of these x a 's are of the 
same order of magnitude. We shall label any one of the latter regions by S, and the former 
regions by R[12 ■■•!'] = {1 > ii > x 2 >■■•> ^ > and its permutations. 

The x-dependence of YliDi in R[12 ■■■£'] is of the form Y\ e a=1 x~ ma , so the x-integral 
encountered in (Bl) is 

r dx & r dxe-! f dxi 1 B 

hx^L T 7 ^'" L x^ ~^ (lns) ' (R5) 

J— Xgl JX£l Xfl_-y J X2 Xi i> 

where M = Yla=i{ m a — 1) and B is determined by how many times the sum J2 a =i( m a — 1) 
reaches zero by varying b from 1 to £'. Clearly B < £', and the only way for B — £' is to 
have all m a = 1, in which case we will call the Ins dependence of T saturated. Otherwise it 
is said to be unsaturated. For the uncut diagrams, x = an d £' = £■ We see in (2.1) that all 
the diagrams except for the 6th order horizontal ladder diagrams C15 to C20 are saturated. 
For the cut diagrams in (6.1) to (6.3), only C16 c is unsaturated, but this diagram does not 
contribute to the final sum because of its subdominance. 

The integral in an S region is like one in an R region with £' reduced. For example, 
suppose Xfj^2 to Xfj are roughly equal but X\ to xp/ (£" = £' — 2) are strongly ordered 
as in R[12 ■■■£"] = R". Then the volume element in the last three variables in spherical 
coordinates is r 2 drdfl 2 , so the integration region is effectively R" but with xe" replaced by 
r. Since £" < £, integrals in S will never lead to saturation. 

Cut diagrams are easier to compute than uncut Feynman diagrams for three reasons. 
First and most trivial, a cut line contains a 5(sz) which makes the corresponding '— ' in- 
tegration simpler to carry out. Second and more importantly, the 5-function demands the 
absence of the '— ' momentum across this line, so the cut line is cut also in the sense of 
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(e) (f) (9) 

Figure 13: Sample cut diagrams to illustrate their computation. 

being an open electric circuit. This generally reduces the number of possible flow diagrams 
and makes the corresponding integral easier to saturate. For example, each of the two-loop 
diagrams in Fig. 13 has only one flow path. Thirdly, the flow pattern often leads to vanishing 
cut diagrams. For example, both diagrams 13(c) and 13(d) are zero from the '+' integration 
around the loop (12345) because the '— ' flows around that loop are all in the same direction. 
This accounts for the equality C17 c = C18 c = in eq. (6.1). 

The cost of this simplicity is the presence of tricky diagrams, which are cut diagrams 
that are apparently logarithmic divergent at large k a+ . This happens whenever there is a 
loop with only one arrow present. Arrows on the top fermion line should not be counted 
for this purpose because their approximate propagators l/{sz + ie) no longer carry any k a+ . 
Figs. 13(a) and 13(b) are examples of such diagrams. This apparent divergence is produced 
by the approximation of replacing the inverse propagators of an upper line by sz + ie, thus 
losing some q + factors needed for convergence. To regulate it we must replace the 5-function 
of the cut line by a smeared 5-function, thereby allowing a small amount of ' — ' momentum 
to flow through, and in the process restoring the lost q + factor. If we do so to 13(a), the k + 
integration is no longer divergent, it will have exactly the same flow path as in 12(a) and 
this is indicated in 13(a) by the dotted arrows, so a pole can be taken at the bottom fermion 
line as shown. The result is given by (Bl) to be (g 4 2/2s)/ 2 • §, the extra \ is there because 
the x integration is bounded by the flow path in 12(a) to be between and 1, so only half 
of the 5(sx) is integrated. When multiplied by the QED vertex and normalization factor 
(2s) 2 / (2m) 2 , one obtains Bl c in (6.1). 

There is another way to compute Fig. 13(a). This is to recognize the fact that the cut 
makes it symmetric in lines 1 and 2, so we may replace the diagram by half the sum of it and 
its crossed diagram. Using the sum rule (3.4), a[12] +a[21] = a[l|2], a cut can be produced at 
the bottom line, and the resulting factor —2ni5(^fsqi + ) makes the qi + integration convergent. 
See Fig. 14. Moreover, the result of this '+' integration is —2ni, exactly the same as if we 
were to do it by residue calculus. The extra factor of \ obtained in the last paragrpah now 
emerges because the double-cut diagram on the right of Fig. 14 is the sum of two diagrams. 

Now we come to the tricky diagram 13(b), which unlike the one-loop case is much more 
difficult to compute by regulating the (^-functions. The reason is that there are now two cut 
lines, the relative magnitude of the small ' — ' flows matters, and that produces once again two 
flow diagrams like 12(b) and 12(c). The poles are no longer at the bottom fermion line and the 
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Figure 14: One-loop tricky diagram and its computation. 



computation is no simpler than the uncut diagrams. We must then compute it by the second 
method. With the double cuts in 13(b), it is symmetric in all the variables. Insisting on 
this symmetry, 13(b) is equal to 1/3! times the sum of 6 diagrams, obtained by permuting the 
bottom gluon lines in all possible ways. From the factorization formula a{l; 2; 3} = a[l|2|3] 
used on the bottom line, the result is equal to a diagram with all its fermion propagators 
cut. So 13(b) is given by (Bl) to be 5f 6 (-27ri) 2 / 3 /(47rs) 2 3! = -g 6 / 3 /24s 2 . Incorporating the 
extra QED factor — (2s) 3 /(2m) 2 for the T matrix, we obtain the contribution of CT5 C to 
M = (2m 2 /g 2 s)T shown in (6.1). 

We conclude this appendix by discussing the remaining expressions in eqs. (6.1)-(6.3). 
We shall use the notation (B2 C ) to denote the spacetime contribution to B2 C , etc. Then 
(B2 C ) = (B2), (C12' c ) = (City = (C12), (C16 c ) = (C16), and (C19 c ) = (C19) can be 
obtained from (2.1). The zero of CT9 C in (6.1) is due to the vanishing of the color factor 
as shown in Fig. 11. Next, we compute (C3' c ) = (C3"). The flow path of this is shown in 
Fig. 13(f), with the dotted arrows indicating the small regulating current which is allowed 
to flow only in the direction shown. The calculation is identical to the uncut diagram ||, 
except that the factor — |ln 2 s is replaced by ~(lns)(— 2ni) (—2ni from the cut, Ins because 
only one uncut line is left on top, and | because only half of the 5-function is integrated). 
The result can then be read off from (2.1) to be 

(C3' c ) = (C3" c ) = -g 2 (3i(\ns)±J 2 I 2 , (B.6) 

which gives rise to the expression in (6.3). 

Finally, we must show that C7' c = and Cl" c = 0. The former can be seen from the flow 
path in Fig. 13(g), where the arrows around the small triangle goes around in a closed loop. 
The latter is so because the scalar diagram in Fig. 7 for Cl" c is symmetrical in lines 2 and 3 
but the triple gluon vertex is antisymmetrical in these two lines. 
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